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Abstract 

- -. - The rate for the production of massive fermions radiated off a pair of massless fermions 

is calculated analytically. Combined with the analytic calculation of the corresponding 
virtual contribution one obtains the order a 2 correction to the cross section which is 



induced by the real and virtual radiation of a pair of massive fermions by massless fermions. 

Approximations valid for large and for small masses are derived and shown to agree with 

<^ ' earlier calculations. 

<N 

> 

^ 1. Introduction 

A large arsenal of approximations and expansions has been developed for QCD calcu- 
li \ lations in those cases where masses are large or small compared to the characteristic energy 



of the problem. There is, however, only a fairly limited number of three loop problems with 
arbitrary nonvanishing masses and energies where answers can be obtained in closed analytical 



43 
Qh! form. The results are first of all useful in their own right; in addition they provide valuable 

Q_i! tests of the approximation methods, strengthening our confidence in these methods in those 



cases where no exact result is available. 

In this paper the imaginary part of the three loop amplitude depicted in Fig.|l| will be 
calculated in closed form. It corresponds to the rate for the production of a pair of massless 
fermions through an external current, with additional radiation of a (real or virtual) pair of 
fermions of mass m. The calculation will be performed in the framework of QED; the transition 
to QCD will be accomplished at the end of the paper. The calculational technique is based on 
the fact that all relevant integrals can be performed in two steps: First one evaluates effectively 
a two loop diagram; the result is subsequently integrated with a fairly simple weight function. 
This technique has been used previously |I|, 0] to calculate individually the contributions of 
the two fermion and the four fermion cuts to the imaginary part. There the limit s ^> m 2 
was considered and only logarithmically enhanced or constant terms were kept. In || the 
same approach was used to evaluate the vertex corrections (including of course the appropriate 
counterterms) for arbitrary m 2 /s. In this paper a similar formula is derived for the four fermion 
cut, completing thus the 0(a 2 ) evaluation of this admittedly rather simple diagram. 

*The complete postscript file of this preprint, including figures, is available via anonymous ftp at 
ttpux2.physik.uni-karlsruhe.de (129.13.102.139) as /ttp94-ll/ttp94-ll.ps 
*e-mail: tt@ttpux2.physik.uni-karlsruhe.de 



Figure 1: Characteristic Feynman diagram describing the production of a pair of massless and 
a (virtual or real) pair of massive fermions. 



The bulk of the paper is formulated for QED as a reference theory, with a defined as 
usual at q 2 = and with fermions of unit charge. The definition of m is irrelevant in the order 
under consideration. At the end of the paper the transition to QCD will be performed and the 
MS definition of the coupling constant employed. 



2. Four fermion contribution 

The four fermion cut through the diagram Fig.[l] corresponds essentially to the part of 
the four fermion production cross section where the light fermion q with mass zero is coupled to 
the external virtual photon or Z boson and the heavy fermion Q with mass m is radiated off the 
light fermion through a virtual photon. For the present problem it is convenient to decompose 
the four fermion phase space {qqQQ) integration into the integration over the three particle 
phase space qq{QQ), where (QQ) denotes a system of fixed invariant mass s', and a (trivial) 
two particle phase space integration for the (QQ) system. This leaves the final integration 
over s'. 

The first step corresponds to the calculation of the rate for the decay of a vector boson 
of mass y/s into a vector boson of mass y/P plus a pair of massless fermions qq. The two 
particle system gives rise to the familiar QQ threshold factor /?(3 — j3 2 )/2. One thus arrives at 
the following integral [Q 
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where u = s'/s and x = m 2 /s. This integral can be solved in a straightforward way numerically, 
and after some effort, also analytically in terms of polylogarithms: 
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where 



a = a/1 + 4x , w = Vl -4x. (3) 

This formula is the main result of this paper. The function g R is evidently zero for m > yfs/2, 
corresponding to x > 1/4. It can be expanded in the limit of small masses, x = m 2 /s <C 1. 
Including terms of order x 3 one obtains 
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The exact result for g R and approximations including successively higher orders in x = m 2 / s 
are displayed in FigJ^a. The curve including terms up to m 4 /s 2 already provides an excellent 
approximation. 



3. Vertex correction 

The vertex correction has been calculated in JJ. Through the two particle cut it leads 
to the following contribution to R: 



with 
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where A = (\/l + 4x — l)/v4x . The function ^ is simply twice the real part of the form factor 
_F 2 defined in in the region x = l/4r > 0. The vertex correction is present for large and 
small m 2 /s as well. The leading term in the heavy mass expansion has also been calculated in 
| | employing a completely different technique 



v s 1 ( m 2 22 N 

S ~— 77 ln — + 77 
m z 45 V s 5 , 



(7) 



As previously discussed in || (see Fig. 2 of 0), the heavy mass expansion provides an excellent 



approximation to the full answer from m 2 ^> s even down to the threshold Am 2 



This 



justifies, for example, the use of eq.(0) for the contribution of virtual top quarks to the Z decay 
rate for the full top mass range. The quality of the expansion is demonstrated in Fig.|3| where 
the full analytic result for g v is compared to the approximation in the range s/m 2 < 4. For 
small x = m 2 /s, on the other hand, one obtains (terms up to 0(x)can be found also in |J) 
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(8) 



In Fig.^b the function g v is displayed in the range from x = 1/4 down to 1/4 x 10~ 2 , to- 
gether with approximations including successively higher orders. Again one concludes that the 
expansion up to m 4 /s 2 provides an excellent approximation in the full mass range for nearly 
massless quarks as well as close to threshold. Upon adding q r and q v and performing the same 
approximations as above, one obtains 
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The quadratic and cubic logarithms compensate as is evident from eqs.(f|) and @. A linear 
logarithm whose origin will be discussed in a moment remains. The sum g v + g R is shown 
in Fig.^a, together with the approximations discussed before. Terms proportional m 2 /s are 
absent in the approximations, whence dotted and dashed-dotted curves coincide. The absence 
of these terms has been observed and dicussed in || in the context of QCD. 



I QCD 

It has become customary to formulate the QCD result in the MS renormalization scheme 
which is better adapted to the limit m — > 0. To perform this transition (still in the abelian 
theory) corrections of order a and a 2 must be redistributed. 

With a defined as fine structure constant in the usual way one finds 
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Figure 2: a) The junction g R describing the production of four fermions in the region < x = 
m 2 /s < 1/4. Solid line: exact result; dashed-dotted line: logarithmic and constant terms only; 
dotted line: including m 2 /s corrections; dashed line: including m 2 /s and m A /s 2 corrections, 
b) Corresponding curves for g v describing virtual corrections. 



Figure 3: The function g v describing virtual corrections in the range s/m 2 < 4 (solid curve) 
and the approximation eq. (Q) (dashed curve). 



The dots indicate other contributions of 0(a 2 ) (for example from two photon exchange or other 
fermion loops) which are irrelevant for the present discussion. The fine structure constant is 
related to the coupling constant in the MS scheme at scale /x 2 through 
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The combination g R + g v + ^ In m 2 /s thus gives the contribution from the (real plus virtual) 
radiation of a massive fermion pair to the total rate or cross section in the MS scheme at scale 
fi 2 = s. It is shown in Fig^b together with the three approximations introduced above. 

The group theoretical transition from QED (U(l)) to QCD (SU(3)) is easily performed 
by multiplying the first order correction by 
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and the second order term by 
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Figure 4: a) The function Q Rj rQ V as described in the text. Solid line: exact result; dashed- dotted 
line: logarithmic and constant terms only; dotted line: including m 2 /s corrections; dashed line: 
including m 2 / s andm A /s 2 corrections, b) The function g R + g v + lnx/4 and its approximations 
describing the same result in the MS scheme. 



The second order correction factor in the small m limit 

2 / R v 1, m 2 \ 2 > , N 11 o / , „ N 13\ 

-^ + / + i ]n 7 J = -C(3)--+x 2 (^-lnx-4C(3) + y j (16) 

i / 8 i 2 56, 16 A/ , 136\ , nl 4 , 

is well behaved. The constant term coincides with the corresponding result given in || (see 
also P|), the quartic mass term with the results of 0. 

To summarize: The rate for the production of massive fermions radiated off a pair of 
massless fermions has been calculated analytically. Combined with the analytic calculation of 
the corresponding virtual contribution one obtains the correction to the total cross section of 
order a 2 which is induced by real and virtual radiation of a pair of massive fermions. The 
transition to QCD and to the MS scheme allows the comparison with approximate formulae 
valid for large or small m 2 /s. Earlier results for a 2 m A /s 2 corrections have been confirmed. 
These leading terms provide an excellent approximation to the full anwer. 
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